Reentrant topological transitions in a quantum wire/superconductor system with 

quasiperiodic lattice modulation 
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We study the condition for a topological superconductor (TS) phase with end Majorana fermions 
to appear when a quasiperiodic lattice modulation is applied to a one-dimensional quantum wire with 
strong spin-orbit interaction situated under a magnetic field and in proximity to a superconductor. 
By density-matrix renormalization group analysis, we find that multiple topological phases with 
Majorana end modes are realized in finite ranges of the filling factor, showing a sequence of reentrant 
transitions as the chemical potential is tuned. The locations of these phases refiect the structure of 
bands in the non-interacting case, which exhibits a distinct self-similar structure. The stability of 
the TS in the presence of an on-site interaction or a harmonic trap potential is also discussed. 

PACS numbers: 71.10.Pm, 03.65.Vf, 85.35. Be, 67.85.-d 



Edge states of topologically nontrivial systems have at- 
tracted attention because they are topologically protected, 
that is, they are stable against weak perturbations that 
do not change the topology of their quantum state. Ma- 
jorana surface states can form at the boundaries or vor- 
tex cores of topological superconductors (TSs)ji"^ and 
there has been a large amount of effort to observe such 
states partly because they are expected to be useful in 
realizing quantum computation. The TS can form in 
a one-dimensional (ID) quantum wire with spin-orbit 
interaction (SOI) that is placed under a Zeeman field 
and close to a bulk superconductor^"— or in cold atoms 
in two-dimensional optical lattices with effective gauge 
fields generated by spatially varying laser fields^ among 
others. 

Here we are interested in the effect of spatial inhomo- 
geneity of the system ^ — i — ' imposed as a modification of 
site energy levels, on the realization of Majorana end 
states. Recently Brouwer and coworkers^'* have stud- 
ied two models, the Dirac equation with random mass 
and a ID spinless superconductor, and obtained the en- 
ergy distributions of the end states. The end modes 
of non-interacting ID systems with quasiperiodic po- 
tentials have also been studied in connection with two- 
dimensional integer quantum Hall systemsi^ and quan- 
tum spin Hall systems.^® 

There is also a growing interest in the role of 
the electron-electron interaction in ID conductors with 
SOI and topological materials. ^^"^^ Stoudenmire and 
coworkers^S showed that while an on-site electron- 
electron interaction U reduces the proximity-induced su- 
perconducting gap on a quantum wire attached to a bulk 
superconductor, the chemical potential range of TS with 
Majorana end states is enlarged by C/ > 0. 

In this work we study the effect of a quasiperiodic site 
level modification of a ID lattice. Such a system has been 
realized in cold-atom systems^^i^i and can also be rele- 
vant when the quantum wire is placed on a bulk super- 
conductor which has the lattice (or superlattice) constant 
incommensurate with that of the quantum wire. We find 



that TS with Majorana end states are observed in sev- 
eral regions with finite widths of the chemical potential 
closer to the band center, which resembles the multichan- 
nel casei ^^i^^'^^ Moreover those TS regions are broadened 
by C/ > 0. This result paves a way to the observation and 
manipulation of Majorana end states in various ID TS 
systems with inhomogeneities, which may be tunable or 
unavoidable. 

Setup. — We study a tight-binding ID fermion model by 
the density-matrix renormalization group (DMRG)i^i^ 
Up to 160 eigenstates of the reduced density matrix p, 
with the sum of discarded eigenvalues of p at each step 
in the last finite-size system loop being typically less than 
10~^, have been retained in the DMRG calculation. 
We adopt the following Hamiltonian, 
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Here, Ca-.i annihilates a fermion with spin a{— at 

site l{— 0,1, . . . , L — 1), hg-.i = c\ iCu^i, t is the nearest- 
neighbor hopping, U is the on-site interaction, A is the 
coupling to the bulk superconductor, a is the Rashba- 
type SOI, Vz is the Zeeman energy, and /i is the chemical 
potential. In the following we set L = 200 and t = 1. We 
take (A, a, Vz) — (0.1, 0.3, 0.3) unless noted otherwise, as 
in Figs. 2, 3 and 6 of Ref. [lO- eo-,; = e; is the site energy 
for spin a on site I. 

For an infinite size system with U = A ^ and 
fo-./ = 0, it is straightforward to obtain the single-particle 
dispersion relation as a function of the quasimomentum 
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E^(k) = t(l - cos(fc)) ± ^Ja^ sm'^{k) + V^ 
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In this Rapid Communication we call them the upper 
and lower Rashba-Zeeman (RZ) bands. 

If the Hamiltonian can be mapped to that of a spinless 
system, the TS state is realized by the introduction of 
the pairing A. When A Vz such mapping is possi- 
ble if /i is in only one of the RZ bands .^°'^^'^° For more 
general discussion on the origin of the topological states 
we refer to Ref. |2^. Note that even when U = 0, while 
only quadratic terms of annihilation and creation opera- 
tors appear in the Hamiltonian, A 7^ introduces non- 
zero matrix elements between states whose number of 
fermions differs by two. The dimension of the Hilbert 
space grows exponentially as the number of lattice sites 
L is increased, strongly limiting the availability of the 
exact diagonalization approach. 

Suppose that we have a lattice system having two Ma- 
jor ana end modes 71,72 such that 7/^ = 71 ± 172 is 
a fermionic operator satisfying {rj'^y — (71 -t- ij2V — 
71 — *72 = ■ If the Majorana operators can be ap- 
proximated by linear combinations of the single-particle 
annihilation and creation operators, 7^ = J2a i(^a\'^<^,^ + 
(''ct ;) ~ li 2), we can think of two single-particle 

wavefunctions \ = J2a li'^alT'^t which |) is the 

empty state, as the "Majorana wavefunctions" . For the 
ground state many-body wave functions in the sectors 
of total number of fermions being even (e) and odd (o), 
|\E'e,o) with energies £^6,0, we can obtain the values of 

„(J = 1,2).20 



(2) 



(*o|cI,,|*e> + (*c|4j*o) 



a 



(*o|<,|* 



(3) 
(4) 



In I20I the phase diagram for ?7 > and U = are ob- 
tained by calculating 1 5*0,0) and checking if the following 
three conditions are met: (i) AE = — Eq vanishes, 
(ii) The left reduced density matrices of the system, ob- 
tained from the density matrices |5'c(o))(^e(o)l by tracing 
out all sites in the right half, have degenerate eigenvalue 
spectrum, (iii) {a]^\} and {a^^]} are spatially localized 
to the different ends. Now we apply these conditions to 
the case with spatial inhomogeneity. 
Quasiperiodic site potential. — We study the effect of a 
quasiperiodic site potential, which is given by 



= Vqcos(k(; - ^c) + <5), 



(5) 



in which Vq > 0, Ic = {L — l)/2 and the phase is (5 = 
unless noted. We choose n — 27rg, in which g = \/5 — 2. 

In the non-interacting case (C/ = A = 0), we can con- 
sider a periodic lattice with gn — Fn-3/Fn, in which 
Fn is the n-th Fibonacci number, to obtain the approx- 
imate eigenenergy distribution. The Fourier transform 
of the site potential then has only components with 
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FIG. 1. (Color online) Differences in energy between the 
ground states in even and odd number sectors of the Hamil- 
tonian plotted against the chemical potential jU. The energy 
spectrum of upper and lower RZ band eigenstates for the 
L = 233 case has been plotted in two lines above the hori- 
zontal AE = line for Vq = 0, and the single particle energy 
spectrum of the A = case for L = 200 has been plotted 
above the AE = line for each value of Vq = 0, 0.2, 0.5, 0.7, 1 
and 1.2. The other parameters are (A,q,V'z) = (0.1,0.3,0.3) 
and U — 0. Insets: (a) The single particle energy spectrum 
of the A = case plotted for k £ [0, 2tt] and Vq = 1, and (b) 
a closeup of the main plot for fi G [1.26 : 1.38] for Vq — 1.2. 



k — ±27r(7„, ±27r(l — gn) but mixes states between the 
upper and lower RZ bands because the spin composition 
of the states in RZ bands depends on the wave number. 
gn rapidly converges to 5 as n is increased, so does the 
eigenvalue distribution. As Vq is increased, the number 
and widths of the gaps in the eigenvalue distribution both 
increase. For Vq ^ t the spectrum exhibits a distinct self- 
similar structure when k is changed, which resembles two 
Hofstadter butterflies^Si^ shifted in energy and overlap- 
ping each other, as shown in the inset (a) of Fig. [TJ We 
note that, for a — Vz — 0, all single body wavefunctions 
are extended for Vq < t and localized for Vq > t in the 
limit of large system (L — > 00) 

In Fig. [T] we plot AE against /i for several values of 
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Vq for L = 200 sites, as well as the eigenenergies of 
the non-interacting case with the quasiperiodic poten- 
tial. The slope of the A£' as a function of /i is often 
close to ±1 outside the TS phase, which corresponds to 
the fact that iVc — is close to ^1 in such regions of /i. 
The two regions with AE = for Vq = correspond to 
the ranges of fi crossing only one of the two RZ bands. 
The RZ bands are gradually mixed and split into several 
mini-bands as Vq is increased also for the open bound- 
ary condition. While the locations and widths of some 
of the regions with vanishing AE resemble those of the 
mini-bands, not all the regions of the chemical potential 
overlapping with one of the mini-bands have AE — 0. 
While the value of AE outside the plateaus at AE de- 
pends on the choice of S, the locations and widths of the 
plateaus almost do not change when 6 is changed. 

We have observed that the other two conditions of 
Ref. Il^, namely (ii) the localized end states and (iii) 
the degeneracy of reduced density matrix eigenvalues, 
are satisfied in the regions in which AE vanishes, and 
they are not when |Ai5| > 10~*. The amplitude distri- 
bution of the wave function of the end Majorana modes 
for (Vq,/x) = (0.2,0.5) is shown in Fig.E^a). The distri- 
butions of ja'-^^p and |a*^^^p are respectively localized to 
the right and left ends of the system. While for Vq > 1, 
the plots of AE{fi, U — 0) exhibit significantly shorter 
plateaus at AE = 0, the two Majorana modes local- 
ized at the two ends are observed inside such plateaus 
as shown for {Vq,h) = (1.2, 1.36) in Fig. [^b). This is 
a nontrivial observation because for |Vq| > 1 we expect 
localized single-particle states for A = 0, and localized 
states would not support global TSs. The reduced den- 
sity matrix eigenvalues are also degenerate, therefore we 
conclude that a TS with end Majorana fermions is re- 
alized in the AE = plateaus. We have also observed 
that the dependence of the regions with Majorana end 
states on the initial phase S is weak in our system, which 
is similar to the non-interacting case,— while the sign of 
AE between these regions depends on the choice of 5. 

While we easily obtain the Bogoliubov quasiparticle 
energies by retaining the value of A and diagonalizing 
the Hamiltonian in the Nambu spinor space ^ii the corre- 
spondence between the degenerate region and the single 
band region in the quasiparticle spectrum has been ob- 
served to be comparable at best to the correspondence 
between the former and the single band region of the 
mini-bands without A. We believe that not only the 
energy spectrum but also the spatial distribution of the 
states is important in the realization of the TS states. 

In the studies of ID topological phases of free fermions 
in bichromatic superlatticesi^ii^ the topological phases 
appear only at special filling factors between the bulk 
bands. Here we emphasize that in our system of the ID 
quantum wire with strong SOI, placed under a magnetic 
field and having a proximity-induced pairing, the topo- 
logical phases are realized in several finite ranges of /i 
that corresponds to finite ranges of the filling factor, as 
observed in Figs. [T] and ^{U = 0), some of them much 



closer to half filling compared to the Vq = case. This 
is the main result of this Rapid Communication. 



(a) {U, Vq, (i) = (0, 0.2, 0.6) □ 
id^ifi mm 




FIG. 2. (Color online) Top: The values of ^^\a.a,i\'^ are 
plotted against the site I for (a) {U,Vq,^) — (0,0.2,0.6), (b) 
{U,Vq,ij,) = (0,1.2,1.36), (c) {U,VQ,fi) = (0,0.2,0.67448), 
(d) (C/,Fq,m) = (0.2,0.2,0.75743), and (e) {U,VQ,fi) = 
(0.5,0.2,0.85472). In (c)-(e), the value of ^ has been cho- 
sen so that A'e = 155.000. Bottom: The eigenvalue spectra of 
the reduced density matrix for the left half of the system in 
states 'I'c and ^'o are plotted for the parameter sets used in 
(c)-(e). 
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FIG. 3. (Color online) The energy difference between ground 
states in the sectors with even and odd number of fermions 
plotted against the number of fermions for U = 0, 0.2, and 0.5. 
The other parameters are (A, a, Vz) = (0.1, 0.3, 0.3) and Vq = 
0.2. Inset: the same plotted against the chemical potential /i. 

Effect of fermion-fermion interaction. — Next we study 
the effect of an on-site fermion-fermion interaction coex- 
isting with the quasiperiodic site potential modulation. 
In Fig.EJthe values of AE{^i, U) for [/ = 0, 0.2 and 0.5 are 
plotted against the number of fermions iVc in the system 
as well as against /x. 

When U is larger, the amount of increase in ^ required 
for adding the same number of fermions at the same fill- 
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FIG. 4. (Color online) AE plotted against fi for (Vq,Vh) = 
(0,0), (0,0.5), (0.2,0), and (0.2,0.5). The other parameters 
are {A,a,V^,U) = (0.1,0.3,0.3,0.5). Inset: AE plotted 
against /j, for Vq = and Vk ~ 0, 0.1, 0.3, 0.5 and 1. 

ing factor becomes larger because of the stronger repul- 
sive interaction. The A£'(/i) = plateaus in the inset 
of Fig. |3] are observed in broader ranges of the chemical 
potential, as expected from Ref. 20 with Vq = 0. Fur- 
thermore, we observe that the plateaus are broader for 
larger U also in terms of the number of fermions^^. 

For a fixed iVc, chosen so that AE approaches zero as U 
is increased, in Fig.[2](c)-(e) we observe that the distribu- 
tions of \a''J\ I become localized toward the ends. We also 
observe that the difference between the eigenstate spectra 
of the reduced density matrices for and , plotted 
in the bottom part of Fig. [2l becomes smaller as AE be- 
comes smaller. Similar behavior of the distributions of 
{a^^ll and the eigenstate spectra of the reduced density 
matrices are observed in other plateaus of AE{fi, U). 
Effect of a trapping potential. — In many cold atom 
experiments, atom clouds are trapped in the vacuum 
by (magneto)optical potentials that are better approx- 
imated by a harmonic or Gaussian potential rather than 
a flat, box-like potential. Also in condensed matter sys- 
tems, the shape of the potential for electrons in the quan- 
tum wire would depend on how it is fabricated. Studying 
the effect of an additional trapping potential on the real- 
ization of the Majorana end fermions is therefore impor- 
tant and has already been conducted for non-interacting 
casesji^i^ To complete this Rapid Communication, we 
study the effect of a harmonic potential, Vu[{l — lc)/lc]'^ 
(I = 0, 1, . . . , L — 1 = 2/c) added to the site potential ([S]). 

Figure m shows AE plotted against fj, for (VqjVh) = 
(0, 0), (0, 0.5), (0.2, 0), and (0.2, 0.5) with (A, a, V^, U) = 
(0.1,0.3,0.3,0.5). The first three plateaus sX AE^Q for 
(Vq, V-h) = (0.2,0) are in the plot. We observe that be- 
tween these plateaus AE plotted for (Vq, Vh) = (0.2, 0.5) 
shows repeated oscillations. 

In the inset of Fig. 0] we have plotted AE for Vq = 



and several values of T^t- F'or increasing T^j the first 
region of vanishing AE' is broadened. The Majorana 
fermion states are localized close to the system boundary, 
where the harmonic trap decreases the effective chemical 
potential, measured from the local site potential, by al- 
most the depth of the trap potential Vh. This explains 
the increase of the value of /i at the upper boundary. 

On the other hand, our system with the quasiperi- 
odic potential has the Majorana state broken down not 
only between the plateau for the system with (Vq, Vh) = 
(0.2,0), but also inside the second and third plateaus, 
when Vh is introduced. We believe that this is because 
the values of the effective chemical potential at which 
the quasiperiodic potential destroys the TS is reached at 
some location even for greater values of /i for Vh > 0. 
This effect seems to compete with the decrease of the 
effective chemical potential at the edge of the fermion 
distribution and AE oscillates. While for simplicity we 
have shown here the results for the case with U = 0.5, the 
discussion above is also consistent with those for other 
values of [/ > 0. 

In conclusion, we have studied the effect of spatial in- 
homogeneity realized by a quasiperiodic site potential 
modulation applied on a tight-binding model of a TS, 
which is a ID conductor with SOI in the proximity of a 
bulk superconductor and under a magnetic field. When 
the modulation is induced, the topological phase appears 
not only when the band is almost empty or almost full, 
but also in several regions of the filling factor (or the 
chemical potential) with finite widths much closer to half 
filling, even when the modulation is strong enough to 
turn the system without SOI and pairing insulating. We 
have also studied the effects of on-site fermion-fermion 
interaction and a harmonic trap. Without the trap, as 
the interaction becomes stronger the TS phase becomes 
wider in the phase diagram, as in the case without the 
spatial inhomogeneity,— in terms of both the chemical 
potential and the number of fermions in the system. 
Our results reveal the possibility of realizing Majorana 
end states in ID systems such as cold-atom systems and 
superconductor-metal heterostructures with inherent or 
imposed inhomogeneties. 
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Recently, we became aware that the interplay of disor- 
der and correlation in ID TSs has also been investigated 
in Ref. [M 
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